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A novel method is presented for the localization of multipoint loss-inducing perturbations in a distributed 
fiber-optic sensor. The proposed simple technique is based on measurement of the transmitted and the 
Rayleigh-backscattered powers of an unmodulated light launched into a sensing fiber. The positions of 
consecutive perturbations are determined by measuring the slopes of the dependence of normalized 
Rayleigh-backscattering power versus the square of normalized transmitted power. It is shown that 
these slopes uniquely depend on the positions of the disturbances along the test fiber. The method 
allows localization of any number of the perturbations that appear one after another at different positions 
along the test fiber without ambiguity. Good agreement is obtained between calculated and experimen- 
tally measured slopes for a loss that was consecutively induced near the source and remote ends of 
2.844-km-long fiber, © 2003 Optical Society of America 
OCIS codes: 060.2370, 290.5870, 060.2430. 



1. Introduction 

The highest state of the art in optical sensing is 
achieved with optical fiber distributed sensors that 
allow measurement of a desired parameter along the 
test fiber. The regions where perturbations occur 
are usually localized by means of optical time-domain 
reflectometry (OTDR) or optical frequency-domain 
reflectometry. 1 - 4 All these methods utilize time- or 
frequency-modulated light sources that allow us to lo- 
calize simultaneously a number of perturbations along 
the test fiber. Meanwhile, for some applications, it is 
important to detect and localize a rare but hazardous 
alarm condition that typically occurs as a single infre- 
quent event, such as a pipe leak, fire, or explosion. 

For such applications recently we proposed a novel 
simple and inexpensive measurement technique 
based on so-called transmission-reflection analysis 
(TRA). The TRA technique utilizes an unmodulated 
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light source, power detectors, and a sensing fiber. A 
novel principle of the localization of loss- induced per- 
turbation, based on measurement of transmitted and 
reflected 5 or Rayleigh-backscattered 6 powers, was 
used in TRA-distributed sensors. Localization of a 
strong disturbance was demonstrated with a maxi- 
mum localization error of a few meters along a single- 
mode sensing fiber a few kilometers long. 

Using the algorithm from earlier research, 5 - 6 we 
can localize only a single perturbation, which may 
appear along the sensing fiber. In practice the num- 
ber of perturbations may affect the sensing fiber ac- 
cidentally or from alarmlike disturbances. Usually, 
however, the probability of alarm conditions being 
detected is low, and therefore they will not appear at 
the same time. 

In this paper we present theoretical and experi- 
mental evidence that the transmission-reflection 
analysis can be modified for detection and localization 
of a number of perturbations that appear one after 
another at different positions along the test fiber. 

In Section 2 we show that localization of a single 
perturbation can be performed with a modified TRA 
method by measuring the slope of dependence of nor- 
malized backscattering power versus the square of 
normalized transmitted power. In Section 3 we 
present mathematical proof that localization of any 
number of consecutive perturbations can be done 
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with a modified TRA method by measuring the cor- 
responding slopes. An experimental verification of 
the proposed technique also is presented. 

2. Single-Perturbati n L calizati n 
First, let us consider modification of the TRA method 
for localization of a single perturbation by using the 
slope of dependence of normalized Rayleigh- 
backscattering power on the square of normalized 
transmitted power. Generally, the TRA method is 
based on the unique relationships between normal- 
ized transmitted and Rayleigh-backscattered powers 
for different locations of the loss-induced disturbance 
along the sensing fiber. 6 Indeed equal losses equally 
decrease the transmitted power for any location of 
disturbance, but the value of Rayleigh backscattered 
power strongly depends on the disturbance location. 

Let us consider the configuration with two plain 
fiber sections separated by a short fiber piece affected 
by a monitored condition (see Fig. 1). Plain fiber 
sections possess Rayleigh scattering and attenuation 
due to light absorption. The power reflection coeffi- 
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Fig. 1. Test fiber configuration for single perturbation: trans- 
mission of the loss-inducing segment; r 1( r. 2t reflections from source 
and remote ends. 



tiple scattering in both directions. The reflections 
with coefficients r l and r 2 from the fiber source and 
remote ends, respectively, must be taken into ac- 
count because even a weak reflection can be com- 
parable with the backscattering. However, we can 
assume that r l 2 « 1 and neglect multiple reflec- 
tions as well. In this case the relationship between 
the normalized transmitted T norm and Rayleigh- 
backscattered fl n orm powers for single perturbation 
can be expressed as 6 



(S„ + r.HR^n - 1) ~ fln 0rm (S Q - r 2 )exp( - 2aL) + S u exp(-2a/ 1 ) 
S u texp(-2aZ!) - exp(-2aL)] + r 2 exp(-2aL) 



(3) 



cient of each Rayleigh-scattering fiber segment can 
be calculated as 7 * 8 

RiL = S(a,/2a)[l - exp(-2a8L)], (1) 



where L is the total sensing fiber length, Z, is the 
distance from the source end to the disturbance loca- 
tion, and S u = S(a,/2a). The normalized transmit- 
ted T norm and backscattered R wmn coefficients in Eq. 
(3) are defined as 



* norm m 1 ' 



T (4) 



T 

•* nr. 

_ J?_ _ S u + r t - (S Q - r 2 )t* exp(-2cxL) - SJ1 - f 1 2 )exp(-2af 1 ) (g) 
"° rm ~~ ^max ~~ S a + r 1 -(S a -r 2 )exp(-2aL) * 



where a s is the attenuation coefficient due to Ray- 
leigh scattering, a is the total attenuation coefficient 
of the test fiber, 8L is the length of the fiber segment, 
and recapture factor S for the fiber is defined as 9 

S = b[n x 2 -ni)ln x \ (2) 

where b depends on the waveguide property of the 
fiber and is usually in the range of 0.21-0.24 for the 
single-mode step-index fiber 9 and n x and n 2 are the 
refractive indices of the fiber core and cladding, re- 
spectively. 

The short fiber piece is affected by monitored con- 
ditions that introduce additional light losses. The 
transmission of a short fiber piece is t y ^ 1. Let us 
assume that the scattering is relatively weak and a 
portion of the scattered light is very small. This 
allows us to simplify the analysis, neglecting mul- 



where T max is the maximum transmittance of the 
initially undisturbed sensing fiber when t x = 1, 

T max = exp(-aL), (6) 

and R max is the maximum backscattering coefficients 
of the undisturbed optical system: 

* m ax - S 0 + r, - (S a - r 2 )exp(-2aL). (7) 

The slope of dependence of normalized backscattering 
power R norm versus the square of normalizedtransmit- 
ted power T norm 2 can be found from Eq. (3.) as 

f ^norm _ 

S tt [exp(-2aii) - exp(-2aD] + r 2 exp(-2aL) 
S„ + r, - (S u - r 2 )exp(-2aL) ' 
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As we can see this slope uniquely depends on per- 
turbation location l v Therefore the location of the 
single perturbation can be found from the experimen- 
tally measured slope as 



all fiber ends was employed to reduce backreflections. 
Standard power detectors were used to measure the 
transmitted and the Rayleigh-backscattered powers. 
To take into account the output laser-power fluctua- 
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The relationship between normalized Rayleigh- 
backscattered power R norm and the square of normal- 
ized transmitted power T nor J is almost linear for a 
single perturbation that affects the test fiber in any 
location [see Eq. (3)]. Figure 2 shows the result of 
the numerical calculation of these relationships when 
additional losses occur at distances l x n = nAL from 
the source end of the test fiber, where n = 0, 1, . . . , 
10, and the interval between bending locations AL = 
284.4 m. Transmitted and backscattered powers 
were normalized with respect to their initial undis- 
turbed values. A typical value for b = 1/4.55 for 
single-mode fiber 10 was used in the calculations. 

The schematic diagram of a TRA-based fiber-optic 
sensor is shown in Fig. 3. Continuous-wave light 
emitted by an amplified spontaneous emission (ASE) 
optical fiber source operating near a 1550-nm wave- 
length with a linewidth of a few nanometers was 
launched into a standard single-mode step-index fi- 
ber 2.844 km long through a 3-dB coupler. The 
launched optical power was -1.1 mW, and the atten- 
uation coefficient of the test fiber a, which was mea- 
sured with an OTDR, was equal to 0. 19 dB /km. The 
optical isolator was used to cancel the influence of 
backreflections on the ASE source. An immersion of 
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Fig. 2. Relations between normalized Rayleigh-backscattered 
power and the square of normalized transmitted power when ad- 
ditional losses occur at distances, l l n =n±L, from the source end 

of the test fiber, where n ~ 0, 1 10 and the interval between 

bending locations, XL = 284.4 m: O, experimental results; 
solid lines, theoretical dependencies. 



tions, a source power meter was utilized. To induce 
the bending losses in the sensing fiber, we used bend- 
ing transducers, also shown schematically in Fig. 3. 
By tuning the bending transducer, we changed the 
normalized transmitted power from its initial undis- 
turbed value equal to 1 to more than -30 dB. The 
bending losses were induced near the source end and 
the near remote end of the test fiber. Good agreement 
between experimental data and theory was obtained 
for (c£ s /a) = 0.68, which is quite reasonable for the fiber 
with total attenuation coefficient a = 0.19 dB/km. 10 
Reflections from the source end and the remote end of 
the sensing fiber, which are equal to 4.7 x 10" 6 and 
1.5 x 10" , respectively, in our experiment, were also 
taken into account in the calculations. 

Experimentally measured slopes for the bending 
losses, which were induced near remote and source 
ends of the test fiber, are equal to 0.109, and 3.63, 
respectively. These values agree well with the val- 
ues calculated with Eq. (8), which are equal to 0.108 
and 3.631, respectively. 

The accuracy of the localization of excess loss with 
the TRA method strongly depends on the value of the 
induced loss. With the TRA method it is easier to 
localize a strong perturbation, but the localization of 
a weak perturbation requires higher accuracy in 
measurements of the transmitted and Rayleigh- 
backscattered powers. In contrast to this the accu- 
racy of the localization of loss with the standard 
OTDR technique depends mainly on the duration of 
the optical test pulse and is practically independent 
of the value of loss. Figure 4 shows the relative 
localization error versus the induced excess loss for 
the TRA method calculated by using the data pre- 
sented in Fig. 2. The relative localization error was 



bending transducer 

Jl. 



ASE source 



test fiber^^^d=5mm 



isolator 



Ji-OQQ 



coupler f 



Reflected power 
meter 



-test fiber 2.844 km - 



Transmitted power 
meter 



Source power 
meter 



Fig. 3. Schematic diagram of the TRA fiber-optic sensor. 
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Fig. 4. Relative localization error versus excess loss for the TRA 
method. 



determined as the localization error divided by their 
value for the strong induced loss that decreased the 
fiber transmission by more than 30 dB. As we can 
see the localization error for weak loss significantly 
exceeds the localization error for strong loss. In the 
experimental measurements in our previous experi- 
ment 11 the localization errors for the strong pertur- 
bation were equal to ±1 m for the source and remote 
ends of a 2.844-km fiber. Therefore localization er- 
rors for 0.1- and 1-dB induced loss are equal to ±22 
and ±3 m, respectively. 

The measurement range and the maximum allow- 
able loss for the TRA and the OTDR methods are 
likely to be the same, because both methods measure 
the backscattering power and the maximum range 
for both methods is namely restricted by the attenu- 
ation of test fiber. 

Summarizing this section, we conclude that the 
modified TRA method allows us to localize the single 
perturbation by measuring the slope of dependence of 
the normalized Rayleigh-backscattered power on the 
square of the normalized transmitted power. 

3. Localization of Multipoint Perturbations 

Let us now verify that any number of consecutive 
perturbations can be localized with the modified TRA 
method. The proof is done by mathematical induc- 
tion. In our analysis we consider the test fiber with 
the same properties and parameters as in the previ- 
ous case. But now we use a configuration with a 
number of plain Rayleigh-scattering fiber sections 
separated by a number of short loss-inducing fiber 
pieces with transmissions t t ^ 1 (see Fig. 5). Let us 
assume that, according to the principle of mathemat- 
ical induction, we have already determined the val- 
ues and locations of the first n perturbations and 
demonstrated that we can find the position of the 
next in -r l)th perturbation without ambiguity. 
Here we consider only the perturbations that appear 
one after another at different positions along the test 
fiber. So, at this moment all initial n perturbations 
induce fixed known losses at known locations and 
only a new (n + l)th perturbation can modify the 
reflectivity and transmittance of the test fiber. 
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Test fiber configuration for multipoint perturbations: 
. transmission of initially disturbed loss-inducing segments; 
* T , transmission of the currently disturbed segment; r,, r 2 , reflec- 
tions from source and remote ends. 
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Because we know the positions and values of all n 
initial perturbations, we can number these according 
to their positions along the test fiber (see Fig. 5) . We 
also suppose without loss of generality that a new 
perturbation is located at a distance l x between the 
kth and (k + Dth initial perturbations, where k is 
unknown. Assuming that the scattering and reflec- 
tions from the fiber ends are relatively weak, we ne- 
glect multiple scattering and reflections in both 
directions. In this case the derivative of the normal- 
ized backscattering power with respect to the square 
of the normalized transmitted power can be ex- 
pressed as (see Appendix A) 



dR n 



S„ I] ;/[exp(-2a/,) 



-exp(-2a/». 1 )] + S 0 2 [«p(-2o/,-) 



j'O 



J~0 



(10) 

where T L = exp(-aL) is the transmission coefficient 
of the test fiber with length L and lj are the location 
of the initial loss-inducing short segments (see Fig. 5). 

This derivative (or the slope of the dependence of 
normalized backscattering power versus the square 
of the normalized transmitted power) depends only 
on the backscattering from the plain segments lo- 
cated on the right side of the unknown perturbation. 
The plain segments located on the left of the un- 
known perturbation do not affect the slope. 

Let us induce for0<£<rc + la helper function 
F(k, n): 



S.S [exp(-2aZi) 



#ma* I! l J 2 '"* 

-exp(-2a/,.0]Il0 2 + ^V 2 n^). 

(11) 

The helper function has a structure similar to the 
expression for the slope [see Eq. (10)], and it de- 
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Fig. 6. Preliminary localization of the {n + l)th perturbation with 
the helper function F(k, n \. 



creases with k for any n (see Fig. 6). The contribu- 
tion in the slope due to backscattering from segment 
Vx* lk+ ii> which is associated with the term 
[exp(-2a/ r ) - exp(-2a/* +1 )] in Eq. (10), is less than 
the possible contribution owing to backscattering 
from the full segment [l ki Z A+l ], which is associated 
with the term [exp(-2a/ A ) - exp(-2a/ A+1 )] in Eq. 
(11). Comparing Eqs. (10) and (11), we conclude 
that, if the measured value of the slope for unknown 
perturbation satisfies 



F(k* + 1, n) < 



3(W) 



<F(k* t n), 



(12) 



the unknown perturbation is located between the 
k*th and (k* + l)th initial perturbations (see Fig. 6). 
Note that, if the measured slope of the unknown per- 
turbation is equal to F(Q, n), the new perturbation 
affects the testing fiber near the source end of the test 
fiber. If the slope is equal to F(n + 1, n) = 
exp(-2aL)r 2 / J R max , the unknown disturbance is lo- 
cated near the remote end of the test fiber. 

Finally, the exact location of the sought-for short 
loss segment can be found: 
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Fig. 7. Relations between normalized Rayleigh-backscattered 
power and the square of the normalized transmitted power for the 
bending losses consequently induced near the remote and source 
ends of the test fiber. 



The increase in the losses leads to a decrease in trans- 
mitted power (line A in Fig. 7). When the square of 
the normalized transmittance decreases to a value 
equal to 0.241 of its initial undisturbed magnitude, 
we stop to increase the bending losses. Afterward, 
keeping constant losses near the remote end, we in- 
duce additional losses near the source end of the test 
fiber. This loading continues until the value of the 
square of normalized transmittance decreases to 
0.061 (line B in Fig. 7). Then, keeping the same 
value of losses near the source end, we gradually 
remove the losses near the remote end of the test fiber 
(line C in Fig. 7). Finally, when the losses near the 
source end are eliminated, all the parameters return 
to their initial undisturbed values (line D). 

All the experimental dependencies presented in 
Fig. 7 possess linear behavior. Experimental data 
show good agreement with theoretical prediction. 
Namely, experimental dependencies A and D practi- 



F(k* + l,n) 



j'k+1 



(13) 



In Section 2 we presented a method for single- 
perturbation localization with a modified TRA tech- 
nique. Now we have demonstrated the algorithm 
for localization of a new perturbation when the values 
and locations of all initial n perturbations are known. 
Therefore, according to the principle of mathematical 
induction, we have demonstrated that a modified 
TRA method can be implemented for the localization 
of any number of consecutive perturbations. 

Figure 7 shows the experimental dependencies of 
normalized Rayleigh-backscattered powers versus 
the square of normalized transmitted powers for the 
bending losses consequently induced near the remote 
and source ends of the test fiber. Measurements 
were performed as follows. Initially the perturba- 
tion occurred near the remote end of the test fiber. 



cally coincide with the calculated dependencies when 
Eq. (3) is used (see Figs. 2 and 7). Experimentally 
measured slopes for lines A, B T C, and D, equal to 
0.109, 3.63, 0.109, and 0.955, respectively, agree with 
values calculated from using Eqs. (8) and (10), equal 
to 0.108, 3.631, 0.108, and 0.957, respectively. Note 
that the slopes were the same for loading and unload- 
ing dependencies. The localization errors that were 
estimated from the difference between the measured 
and the calculated slopes do not exceed 2 m for any 
location of perturbation. Note, however, that in 
practice different noise origins and system imperfec- 
tions such as temporal drifts of fiber and photodetec- 
tors parameters, additional uncontrolled losses, etc. 
may also contribute to a decrease in accuracy. 
The TRA method can be implemented for localiza- 
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Fig. 8. Relations between normalized Rayleigh-backscattered 
power and the square of the normalized transmitted power for two 
perturbations, A + D. synchronously and, A and D, independently 
induced near the remote and source ends of the test fiber. 



tion of any number of consecutive perturbations that 
occur one after another along the test fiber during the 
monitoring period. In contrast to the OTDR the pro- 
posed method cannot be used for localization of the 
perturbations in a fiber optical system already in- 
stalled and that already has many induced losses. 
The other natural question to be addressed about the 
TRA method is the following: What happens when 
two perturbations affect the test fiber simulta- 
neously? To answer this question, consider calculat- 
ing the dependence of normalized Rayleigh- 
backscattered power versus the square of the 
normalized transmitted power for two equal pertur- 
bations that induce losses near the source and remote 
ends at the same time (curve A + D in Fig. 8). The 
dependence exhibits clear nonlinear behavior. As 
shown above for any number of consecutive pertur- 
bations this dependence should be linear. Figure 8 
also shows the normalized Rayleigh-backscattered 
power versus the square of the normalized transmit- 
ted power for perturbations that affect the testing 
fiber one after another near the remote (line A) and 
source (line D) ends. Both of the latter dependencies 
exhibit clear linear behavior. The nonlinear behav- 
ior of the dependencies of normalized Rayleigh back- 
scattered power versus the square of normalized 
transmitted power indicates that the testing fiber is 
affected by two or more perturbations simulta- 
neously. When the particular dependence shown in 
Fig. 8 (curve A + D) is used, it is possible to localize 
at least one perturbation. Indeed the value of the 
normalized Rayleigh backscattered power at the 
point at which normalized transmitted power is equal 
to zero directly shows the location nearest to the 
source-end perturbation. Therefore, analyzing the 
curve A + D, we conclude that two perturbations 
affect the test fiber simultaneously and one of the 
perturbations induces losses near the source end. 
Nevertheless a complete analysis of the different pos- 
sibilities, even for two synchronous events, is notice- 
ably complicated and beyond the scope of this paper. 



4. Summary 

We have demonstrated that the modified TRA 
method can be implemented for the localization of 
any number of loss-inducing perturbations that ap- 
pear one after another at different positions along the 
test fiber. The localization of the perturbations with 
a modified TRA method was performed by measuring 
the slope of dependence of the normalized Rayleigh- 
backscattered power versus the square of normalized 
transmitted power. It was shown that these slopes 
uniquely depend on the position of the perturbations 
along the test fiber. We believe that the proposed 
TRA technique should be attractive in the eventual 
realization of a compact and inexpensive distributed 
fiber-optic sensor. 

Appendix A: Calculation of the Slope 

Let us consider an optical system that consists of a 
number of Rayleigh-scattering plain segments sepa- 
rated by a number of short loss-inducing fiber pieces 
(see Fig. 5). The transmittance of initial n loss- 
inducing short segments located at distances lj from 
the source end is L s 1, where j = 1, n. The trans- 
mittance of an unknown loss-inducing segment that 
is located at distance l x is t x . Neglecting multiple 
scattering and reflections in both directions, we can 
calculate the normalized power-reflection coefficient 
of the optical system as 

R } + Rn + i? 3 
-»*norni p > (Al) 

max 

where R l corresponds to Rayleigh backscattering 
from segments placed on the left of the unknown 
perturbation plus the Fresnel reflection from the 
source end of the test fiber: 

Ri = r x + S u 2 \T t *{\ - exp[-2a(/, + 1 

R 2 is the Rayleigh-backscattering coefficient of two 
segments that are placed around the unknown per- 
turbation: 

R 2 = S u 11 - exp[-2a(/ x - /,)]} 

+ T lx 2 {l - exp[-2a(/ w - l x )]}t x \ (A3) 

and R 3 is the power-reflection coefficient for segments 
that are placed on the right of the disturbed region 
plus the term related to the reflection from the re- 
mote end of the test fiber: 

Rs = t x 2 ls a 2 r /t 2 {l-exp[-2a(Z l+1 -Z ( )]} 

* II '/ + T L 2 r2 f[ </ > (A4) 

y-0 j~Q / 
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where T H - exp(-a/,) is the transmission coefficient of 
the fiber segment with length (see Fig. 5) and r lf r 2 
are the reflections coefficients from the fiber source 
and remote ends, respectively. In Eqs. (A2)-(A4) we 
assign l 0 = 0, t 0 = i f and [ n+% = L. We emphasize 
that the condition * 0 = 1 does not mean that pertur- 
bation cannot appear near the source end of the test 
fiber. The disturbance can affect the fiber near the 
source end, but it should be marked as the first per- 
turbation with transmittance t x at the Zj = 0 distance. 

The normalized transmitted coefficient of the opti- 
cal system, which is affected by all n + 1 perturba- 
tions, is defined as 



T 2 = t 2 TT t 2 

A norm \. J. J * 

j=l 



(AS) 



Note that changes in the first n perturbations are 
almost finished before a new one is started, so that 
tj - constant for 7 = 1, n. Therefore the total nor- 
malized transmitted power for all n + 1 perturba- 
tions T norm 2 can change only because of the change in 
current perturbation. The differential d(T norm 2 ) is 

n 

t(T n0T J) = 3(f, 2 ) II 'A (A6) 

>=i 

The derivative of the normalized backscattering 
power with respect to the square of the normalized 
transmitted power is expressed as 



4- ,_ „ -f 3 



aiT^*) d{T 

norm 



\S„ J! t/[exp(-2al x ) 

r TT t 2 1 j '° 



j=0 



-exp(-2a/* + 1 )] + S« 2 [exp(-2a/ i ) 
-exp(-2*l tn )]Ut* + T L 2 r 2 fltj 



j=0 



(A7) 
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